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Abstract 

The resistance between arbitrary two nodes in a resistor network is 
obtained in terms of the eigenvalues and eigenfunctions of the Lapla- 
cian matrix associated with the network. ExpUcit formulas for two- 
point resistances are deduced for regular lattices in one, two, and 
three dimensions under various boundary conditions including that of 
a Mobius strip and a Klein bottle. The emphasis is on lattices of finite 
sizes. We also deduce summation and product identities which can be 
used to analyze large-size expansions of two-and-higher dimensional 
lattices. 
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1 Introduction 



A classic problem in electric circuit theory studied by numerous authors over 
many years is the computation of the resistance between two nodes in a 
resistor network (for a list of relevant references up to 2000 see, e.g., P). 
Besides being a central problem in electric circuit theory, the computation of 
resistances is also relevant to a wide range of problems ranging from random 
walks (see [2, 3^ and discussions below), the theory of harmonic functions j4j, 
first-passage processes [S], to lattice Green's functions [Hj. The connection 
with these problems originates from the fact that electrical potentials on a 
grid are governed by the same difference equations as those occurring in the 
other problems. For this reason, the resistance problem is often studied from 
the point of view of solving the difference equations, which is most conve- 
niently carried out for infinite networks. Very little attention has been paid 
to finite networks, even though the latter are the ones occurring in real life. 
In this paper we take up this problem and present a general formulation for 
finite networks. Particularly, known results for infinite networks are recov- 
ered by taking the infinite-size limit. In later papers we plan to study effects 
of lattice defects and carry out finite-size analyses. 

The study of electric networks was formulated by Kirchhoff |2| more than 
150 years ago as an instance of a linear analysis. Here, we analyze the prob- 
lem along the same line of approach. Our starting point is the consideration 
of the Laplacian matrix associated with a network, which is a matrix whose 
entries are the conductances connecting pairs of nodes. Just as in graph 
theory that everything about a graph is described by its adjacency matrix 
(whose elements is 1 if two vertices are connected and otherwise), we expect 
everything about an electric network is described by its Laplacian. Indeed, 
in Section 2 below we shall derive an expression of the two-point resistance 
between arbitrary two nodes in terms of the eigenvalues and eigenvectors 
of the Laplacian matrix jSj. In ensuing sections we apply our formulation 
to networks of one- dimensional and two-dimensional lattices under various 
boundary conditions including those embedded on a Mobius strip and a Klein 
bottle, and lattices in higher spatial dimensions. We also deduce summation 
and product identities which can be used to reduce the computational la- 
bor as well as analyze large-size expansions in two-and-higher dimensions. 
In subsequent papers we shall consider large-size expansions and effects of 
defects in finite networks, the latter a problem that has been studied in the 
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past only for infinite networks jU]. 

Let Lbe a resistor network consisting of M nodes numbered hy i = 
l,2,---,7Vr. Let 

^ij — ji be the resistance of the resistor connecting nodes i 
and j. Hence, the conductance is 



so that Cij = (as in an adjacency matrix) if there is no resistor connecting 
i and j. 

Denote the electric potential at the ith node by Vi and the net current 
flowing into the network at the zth node by /j (which is zero if the ith node 
is not connected to the external world). Since there exist no sinks or sources 
of current including the external world, we have the constraint 



Eh 

1=1 



0. 



The Kirchhoff law states 



AT 



1,2, 



where the prime denotes the omission of the term j = i. Explicitly, 

LV = 1 



where 



Cl 
-Cl2 



-Ci2 
C2 



-CiM \ 

-C2N 



J 



and 



Ci 



(1) 



(2) 

reads 
(3) 



(4) 



(5) 



r 

Here V and / are A/ -vectors whose components are Vi and /j respectively. 
We note in passing that if all nonzero resistances are equal to 1, then the 
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off-diagonal elements of the matrix — L are precisely those of the adjacency 
matrix of L. 



The Laplacian matrix L of the network C is also known as the Kirchhoff 
matrix, or simply the tree matrix; the latter name derived from the fact 
that all cofactors of L are equal, and equal to the spanning tree generating 
function for £ [TU] . 

To compute the resistance Ro,^ between two nodes a and /3, we connect 
a and /? to the two terminals of an external battery and measure the current 
/ going through the battery while no other nodes are connected to external 
sources. Let the potentials at the two nodes be, respectively, and Vp. 
Then, the desired resistance is the ratio 

R., = (6) 

The computation of the two-point resistance i?o/3 is now reduced to solving 
(jni) for Va and Vg with the current given by 

/, = /(5,,-M- (7) 

A probabilistic interpretation: 

The two-point resistance has a probabilistic interpretation. Consider a 
random walker walking on the network L with the probability 

Pi^j = Cij/ci (8) 

of hoping from node i to node j, where Pi^j can be different from Pj^i. Let 
P(a; (3) be the probability that the walker starting from node a will reach 
node j3 before returning to a, which is the probability of first passage. Then 
one has the relation [21 E] 

P{a-,P) = -^ (9) 

where Cq, is defined in (0). If all resistances are 1, then Q becomes 

P{a-,P) = ^—. (10) 

where is the coordination number of, or the number of nodes connected 
to, the node a. 
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2 The two-point resistance: A theorem 

Let and k be the eigenvectors and eigenvalues of L, namely, 

L^, = /i^i, i = 1,2,---,AA. 

Let ipia, a = 1, 2, ■ ■ ■ , A/" be the components of \E'i. Since L is Hermitian, the 
\E'i's can be taken to be orthonormal satisfying 

a 

Now the sum of all columns (or rows) of L is identically zero, so one of 
the eigenvalues of L is identically zero. It is readily verified that the zero 
eigenvalue /i = has the eigenvector 



V^i^ = l/VAA, a = 1,2,---,A^. 

We now state our main result theorem. 
Theorem: 

Consider a resistance network whose Laplacian has nonzero eigenvalues 
li with corresponding eigenvectors \E'j = (V'li, V'125 ■ ■ ■ , "^ia/")? i = 2, 3, ■ ■ ■ , A/". 
Then the resistance between nodes a and (3 is given by 



1 



i=2 "« 



2 



= H rki« ~ • (11) 

67 



Proof. 

We proceed by solving Eq. (0) by introducing the inverse of the Laplacian 
L, or the Green's function f^. However, since one of the eigenvalues of L is 
zero, the inverse of L must be considered with care. 

We add a small term e I to the Laplacian, where I is the JV x JV identity 
matrix, and set e = at the end. The modified Laplacian 

L(e) = L + el 

is of the same form of L except that the diagonal elements q are replaced by 
Cj + e. It is clear that L(e) has eigenvalues k + e and is diagonalized by the 
same unitary transformation which diagonalizes L. 
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The inverse of L(e) is now well-defined and we write 

G(e) = L-i(e). 

Rewrite Q as L(e)\?(e) = / and multiply from the left by G(e) to obtain 
V{e) = G(e)/. Explicitly, this reads 

V,{e) = J2Gij{e)I„ t = l,2,---,Af, (12) 

i=i 

where Gij{e) is the ij-th elements of the matrix G(e). 
We now compute the Green's function Gij{e). 

Let U be the unitary matrix which diagonalizes L(e) and L, namely, 

U^LU = A 
UtL(e)U = A(e). (13) 

It is readily verified that elements of U are Uy = tpji, and A and A(e) are, 
respectively, diagonal matrices with elements /j 6ij and (Zj + e) 6ij. 

The inverse of the second line of (jl3|) is 

UtG(e)U = A-i(e) 
where A~^(e) has elements {k + e)~^6ij. It follows that we have 

G(e) =UA-i(e)U^ 



or, explicitly. 



where 



^ + ^7a/3(e) (14) 



9ap{e) = 2^ y— 7- (15) 
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The substitution of (jl4|) into (|12p now yields, after making use of the con- 
straint (HI), 



It is now safe to take the e ^ hmit to obtain 



(16) 



Finally, by combining ©, ((ZI) with (fT^ we obtain 

Rap = 9aaiO) + gfS/siO) - ga/siO) " fi'/3a(0) 

which becomes (jlip after introducing p5|) . QED • 

The usefulness of (II ip is illustrated by the following examples. 

Example 1: Consider the 4- node network shown in Fig. 1 with the 
Laplacian 



/ 


2ci 


-Cl 





-Cl \ 




-Cl 


2Ci + C2 


-Cl 


-C2 







-Cl 


2ci 


-Cl 


V 


-Cl 


-C2 


-Cl 


2Ci + C2 / 



where ci = 1/ri, C2 = l/r2. The nonzero eigenvalues of L and their orthonor- 
mal eigenvectors are 



2 = -(1,-1, 1,-1] 



vl>3 = -^(-1,0,1,0) 



1 

V2' 



vl>4 = — (0,-1,0,1). 



h = 4ci, 
^3 = 2ci, 

U = 2(ci + C2 

Using (fTTj) . we obtain 

1 / \2 1 / n2 1 / 

i?13 = r(^21-V'23 +r(^31-^33 +^(^/'41-^43 
f2 '3 ^ <4 ^ 

1 2 1 2 1 

- (^^21 - V'22) + Y (^31 - ^^32) + — (V^41 - ^^42) 
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n, 

ri(3ri + 2r2) 
4(ri + r2) 



Example 2: We consider AA(A/' — l)/2 resistors of equal resistance r on 
a complete graph of M nodes. A complete graph is a network in which every 
node is connected to every other node. The Laplacian is therefore 



(M-l -1 
-1 M- 1 



complete graph 



\ 



-1 
-1 



-1 
-1 



-1 
-1 



-1 \ 

-I 



U-l -1 
-1 U -I 



(17) 



It is readily verified that the Laplacian (fT7|) has eigenvalues Aq = and 



r ,n 



components 

'4'nc 



1,2,..., A/" — 1, with corresponding eigenvectors having 



exp I 



27ina\ 



n,a = 0, 1, ...,Af - 1. 



It follows from (pT| the resistance between any two nodes a and f3 is 



R. 



■a,l3 



r ^ 1 



n=l 



(18) 



In subsequent sections we consider applications of (fTT|) to regular lattices. 



3 One-dimensional lattice 

It is instructive to first consider the one- dimensional case of a linear array of 
resistors. We consider free and periodic boundary conditions separately. 

Free boundary condition: 

Consider — 1 resistors of resistance r each connected in series forming 
a chain of A^ nodes numbered from 0, 1, 2, ■ ■ -, to A^ — 1 as shown in Fig. 2, 
where each of the two end nodes connects to only one interior node. This is 
the Neumann (or the free) boundary condition. The Laplacian ^ assumes 
the form 

T free —1 rp free 

-^{Afxl} — ' 
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where T^*^*^ is the N x N matrix 



/ 1 -1 
-1 2 -1 



-1 free 
■N 



\ 













\ 




-1 2 -1 
-1 1 



(19) 



The eigenvalues and eigenvectors of T^r can be readily computed (see, for 
example, jH]), and are 



A. 



2(1 — cos (pn 



(N) 
nx 



N 



n = 0, 1,---,A^ 



n = 0, all X 



n = 1,2, 



N-1, all X, (20) 



where 0„ = mr/N. Thus, using (fTT|) . the resistance between nodes xi and X2 



IS 



where 



r [cos {Xi + |)0„ - cos {X2 + l)(f)r. 

N l-cos0„ 
r[F7v(xi + X2 + 1) + Fn{xi - X2) 

-^F^(2xi + l)-^F^(2x2 + l) 



1 ^1 1 - cos 



n=l 



COS( 



(21) 



(22) 



Note that without loss of generality we can take > i < 2N. 

The function Fj~^{€) can be evaluated by taking the limit of / — >■ of the 
function /i(0) — /i(^) evaluated in (jHH) below. It is however instructive to 
evaluate F^{£) directly. To do this we consider the real part of the summation 



1 ^1 1 - e' 



oi<t>n 



First, writing out the real part of the summand, we obtain 



UeTNitj 



1 



N-l 



N ^ 



n=l 



1 — COS (pn — COS i(f)n + COs(£ — 1)0„ 

2(1 - cos(/)„,) 



^[f^(i) + f^(£)-F;v(^-i; 



(23) 



Second, we expand the summand to obtain 

1 N-l l-l 
n=l £'=0 

and carry out the summation over n. The term i' = yields F]^{1) and 
terms i' > 1 can be explicitly summed, leading to 



1 - i-iy 

I _ gini'/N 



< 2N. 



(24) 



We now evaluate the real part of T^i^t) giving by (f^ . Using the identity 

1 ^ 



7^e 



1 - e" 



1 
2' 



< < 27r, 



(25) 



we find 



7^eT^(£) = F^.(l)- — 



2£-3 



-11 



(26) 



Equating (PHj) with and noting FAr(l) = 1 — we are led to the 

recursion relation 



2i-l 



which can be solved (Cf. section 9 below) to yield 

I.I 1 f^^+\i\ 



N 



(27) 



where [x] denotes the integral part of x. The substitution of (j2ZI) into (j2H) 
now gives the answer 



(28) 
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which is the expected expression. 

Periodic boundary conditions: 

Consider next periodic boundary conditions for which nodes and — 1 
are also connected as shown in Fig. 3. The Laplacian ^ of the lattice is 
therefore L|^^^| = r^^ T^^"^ where 

-1 \ 




/ 2 -1 
-1 2 -1 



rpper 
-"-AT 




\ -1 



-1 2 -1 

-1 2 y 



(29) 



The eigenvalues and eigenfunctions of T^"^ are well-known and are, respec- 
tively, 



2 (^1 — cos2(/)j 
1 



N 



n,x = 0,1, - ■ ■ ,N - 1. 



(30) 



Substituting pUll into (fTT|) . we obtain the resistance between nodes Xi and 



X2 



where 



N-l L«2xi<ji„ ^i2x24>r. 



N 2(1 -cos 20, 

rGNixi - X2) 



N-l 



N 



^-1 

n=l ^ 



1 — cos(2£0. 



cos 20, 



(31) 



The function Gn{^) is again evaluated in a special case of the identity 
()62|) below. Following analyses in section 9, one obtains the recursion relation 

which can be solved to yield 



gm = \^\-^Vn. 



(32) 
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It follows that we have 



-R{^xi}(^i5^2) = r |xi - X2I 



N 



(33) 



The expression (j33|l is the expected resistance of two resistors \xi — X2I r and 
[N — \xi — X2I) r connected in parallel as in a ring. 



4 Two-dimensional network: Free boundaries 



Consider a rectangular network of resistors connected in an array of M x 
nodes forming a network with free boundaries as shown in Fig. 4. 

Number the nodes by coordinates {m, n},0<m<M — l,0<n<A^ — 1 
and denote the resistances along the two principal directions by r and s. The 
Laplacian is therefore 



free 



M 



1 free 
-N 



(34) 



where ® denotes direct matrix products and Tj^®^ is given by ()19|) . The 
Laplacian can be diagonalized in the two subspaces separately, yielding eigen- 
values 



and eigenvectors 



2r"^(l - cos6'm) + 2s ^(l-cos0n) 



(35) 



It then follows from (fTT|) that the resistance -Rfree between two nodes ri = 
{xi,yi) and r2 = {x2,y2) is 



,(M),/,(7V) 



M-1 N-1 



-R{MxAr}(ri, r2) - 



^free 



{m,n);{xi ,yi) 



_ J. free 



m=0 n=0(m,n)7^(0,0) 



A 



{m,n) 



r 

'~ N 

M-l N 



Xi - X2 



S 
M 



yi - 2/2 



+ 



MN 



m=l n=l 



1 [cos (xi + l^6m COS + l^4>n- COS (^X2 + l^Om COS (jj2 + \ 



r~l(l — COS ^m) + S ^(1 — COS0„) 



(37) 
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where 

= <Pn = ^• 



Here, use has been made of ()28|1 for summing over the m = and n = terms. 
The resuhing expression (jTTj) now depends on the coordinates xi,yi,X2 and 
y2 individually. 

The usefulness of fjTfjl is best illustrated by applications. Several examples 
are now given. 

Example 3: For M = 5, N = 4,r = s, the resistance between nodes 
{0, 0, } and {3, 3} is computed from (jHTj) as 

«{S;.,({0,0.U3,3}) ^ (i.i.i^). 

= (1.70786...) r. (38) 

Example 4: For M = = 4, we find 

ofree (rr. ni QiA _ (r + g) (r^ + 5rg + s'){3r^ + 7rg + 3s^) 
%x4}(|0, 0|, |3, 3|) - 2(2r2 + 4r. + .2)(r2 + 4r. + 2.2) " ^^^^ 

Example 5: We evaluate the resistance between two nodes in the interior 
of a large lattice. Consider, for definiteness, both M, N = odd (for other 
parities the result below is the same) and compute the resistance between 
two nodes in the center region, 

, /M-1 N-1 ^ 

ri = {xi,yi) = [ — - — +pi, ^ +qi) 

. /M-1 N-1 N 

r2 = {X2,y2> = [ — ^ — +P2,— — + q2), 

where Pi,qi << M, N are integers. The numerator of the summand in ()37|1 
becomes 

[— + Pic'mj cos + - cos +P2fc'mj cos + g20nj 

= cospi6'm cosgi0„ — cosp26'm cosg20n) 5 m, 71 = even 

= (sinpi6'm singi0„ - sinp26'm sing20n) , m,n = odd 

= i^smpiOm cosgi0„ — sinp26'm cosg20n) , ^ = odd, n = even 

= (^cospiOm singi0„ — cosp2^m sing20n) , Tn = even, n = odd 

13 



cos 



and the summation in (|37p breaks into four parts. In the M, N oo hmit 
the summations can be replaced by integrals. After some reduction we arrive 
at the expression 

I r2n r2TT / 1 — COS (xi — X2)6' COs(?/i — ?/2)0' 



f2TT 









{2tiY Jo ^ Jo V r-i(l - cos^) + - COS0) y 

(40) 

Our expression pn|) agrees with the known expression obtained previously 
P]. It can be verified that the expression (pn|) holds between any two nodes 
in the lattice, provided that the two nodes are far from the boundaries. 



5 Two-dimensional network: periodic bound- 
ary conditions 

We next consider an M x network with periodic boundary conditions. The 
Laplacian in this case is 



per 



Irppcr 



{MxN} - ' -LM w Iat + S ^Im ® T^'^'' (41) 



where T^'^ is given by (j^ . The Laplacian PT|) can again be diagonalized 
in the two subspace separately, yielding eigenvalues and eigenvectors 



'4^{m,n)\{x,y) 



2r~^(l - cos2^„) + 2s"" (1 - cos20„) 
1 



^i2x0m ^i2y<f> n 



Vmn 

This leads to the resistance between nodes ri = (xi,?/i) and r2 = {x2,y2 

2 

'4'{m,n);{xi,yi) ~ ^{m,n);(x2,J/2) 



(42) 



pPGr / 



M-1 N-1 

E E 

m=0 n=0{m,n)7^(0,0) 



A 



(m.n) 



r 

N 



[Xi - X2) 

M 



+ 



M 



1/1-1/2 



(2/1 - 2/2)^ 



N 



1 ^ ^l-cos[2{xi-X2)e.^ + 2{y,-y2) 
^ MiV^,^, r-i(l-cos2^„) + s-i(l-cos2<j 



(43) 
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where the two terms in the second hne are given by (j33|) . It is clear that the 
result depends only on the differences \xi — X2I and \yi — 1/2 1, as it should 
under periodic boundary conditions. 

Example 6: Using (PHj) the resistance between nodes {0, 0} and {3, 3} 
on a 5 X 4 periodic lattice with r = s is 

«iT«,({o.o}.{3,3}) ^ {1^1^^!^} 

= (0.680937...) r. (44) 

This is to be compared to the value 1.707863... r for free boundary conditions 
given in Example 3. It can also be verified that the resistance between nodes 
({0, 0}) and ({2, 1}) is also given by (jH)) as it must for a periodic lattice. 

In the limit of M, N ^ 00 with Iri — r2| finite, (1431) becomes 



1 .2^ /■27r 1 - cos [(Xi - X2)^ + (2/1 - 1/2)0 

(/tt)^ Jo Jo r ^(1 — cost^)+s ^(1 — COS0) 
1 f"^"" ^j. f'^'" ^Q ^~ cos(xi - X2)0 cosjyi - ^2)0 ^^^^ 



(27r)2 7o Jo rr^(l - cos6') + s-i(l - cos 
which agrees with (j40|) . 



6 Cylindrical boundary conditions 

Consider an M x resistor network embedded on a cylinder with periodic 
boundary in the direction of M and free boundaries in the direction of N. 

The Laplacian is 

LfMxTV} = T^M ® liv + Im ® T^^^ 

which can again be diagonalized in the two subspaces separately. This gives 
the eigenvalues and eigenvectors 



\m,n) = 2r ^(1 - COs26'rri) + 2S ^(1-COS0„) 
1 

M 



-,/,cyi - pi2xe^ / {N) 
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It follows that the resistance -Rfree between nodes ri = {xi,yi) and r2 
ix2,y2) is 



r 

N 



Xi - X2 



A/-1 Af-1 

E E 



A/-1 Af-1 

E E 

m=0 r!,=0(m,n)^(0,0) 



'^{m,n);{xi,yi) '^{m,n);{x2,y2) 



M 



+ 



M 



yi - y2 



' Ci^ + C2^ ~ 2CiC2Cos2{xi - X2)9„ 



MN ^^^^\ r-i(l-cos2e^) + s-i(l-cos. 



where 



Ci = cos [yi + 



C2 = cos {y2 + -^)(pn- 



(46) 



It can be verified that in the M, N 00 limit (j46p leads to the same expres- 
sion (j4Uj) for two interior nodes in an infinite lattice. 

Example 7: The resistance between nodes {0, 0} and {3, 3} on a 5 x 4 
cylindrical lattice with r = s is computed to be 



^Sx4}({0,0},{3,3} 



3 3 5023\ 

— \ \ h 

10 5 8835 / 

1.46853...) r. 



(47) 



This is compared to the values of (1.70786...) r for free boundary conditions 
and (0.680937...) r for periodic boundary conditions. 



7 Mobius strip 

We next consider an M x resistor lattice embedded on a Mobius strip of 
width and length M, which is a rectangular strip connected at two ends 
after a 180° twist of one of the two ends of the strip. The schematic figure of 
a Mobius strip is shown in Fig. 5(a). The Laplacian for this lattice assumes 
the form 



T Mob _ -1 
■^{MxAT} — ' 



H 



M 



I 



N 



K 



M 



'N 



M 



1 free 
■N 



(48) 
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where 



/ 




K 



N 



H 



N 





VI 



T%'' + K 



1 \ 






y 



/ 




1 \ 

1 





y 



N 



1 

VI 

\ 


2 -1 
-1 2 



are N x N matrices. Now Tj^*'® and J^r commute so they can be replaced 
by their respective eigenvalues 2(1 — cos^n) and (— 1)" and we need only to 
diagonalize an M x M matrix. This leads to the following eigenvalues and 
eigenvectors of the Laplacian (jlH)) [TH IT^ : 



/ 


2 


-1 • ■ 




-1 


2 • ■ 







■ ■ 


V 





■■ 



A 



{m,n) 



,Mob 

T (m,n)\(x,y) 



2r-^cos [(4m + 1 - (-1)" 
1 



TT 



+ 2s"Ml - cos 



exp 



2M 

xn 



niT 
iV 



2M 



(TV) 
ny 



(49) 



where ipl^^ is given in ^lU^ . Substituting these expressions into (fTTj) and after 



a little reduction, we obtain 



oMob 
It 



M-l N-l 

E E 

m=0 n=0{m,r!.)7^(0,0) 



„/,Mob _ „;.Mob 

^{m,n);(xi,y-i) ^ {m,n);{x2,y2) 



A 



(m.n) 



r 
N 



Xi - X2 



[Xi - X2) 

M 



M- 



MN 



E E 



Ci 2 + C2 2 - 2C1C2 COS 


(xi-X2)(4m + l-(-ir)2f^ 




1 - cos (4m + 1 - (-l)")2ij 


+ s ^{1 - COS0„) 



(50) 



where Ci and C2 have been given in 
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Example 8: The 2x2 Mobius strip is a complete graph of = 4 nodes. 
For r = s the expression (|5(J|) gives a resistance r/2 between any two nodes 
which agrees with (fTHj) . 

Example 9: The resistance between nodes (0, 0) and (3, 3) on a 5 x 4 
Mobius strip with r = s is computed from (j5Up as 

= (0.896367...) r. (51) 

This is to be compared to the corresponding values for the same network 
under other boundary conditions in Examples 3, 6, and 7. 



8 Klein bottle 



A Klein bottle is a Mobius strip with a periodic boundary condition imposed 
in the other direction. We consider an M x A resistor grid embedded on a 
Klein bottle, a schematic figure of which is shown in Fig. 5(b). 

Let the network have a twisted boundary condition in the direction of the 
length M and a periodic boundary condition in the direction of the width N. 
Then, in analogous to (j48|) . the Laplacian of the network assumes the form 

Km ® Jivl + s'Hm ® T^"- (52) 



J Klein 
^{MxN} 



H 



M 



I 



N 



Now the matrices Jn and T^'^'^ commute so they can be replaced by their 
respective eigenvalues ±1 and 2(1 — cos20ji) in and one needs only to 
diagonalize an M x M matrix. This leads to the following eigenvalues and 
eigenvectors for L^]^™^| (TH [T2] : 



>(m,n) [ 



I Klein 

T(m,n);(x,y) 



2r^ 



1 - cos ((2m + ^) ^) + 2s"^ (l - cos 



2n'K\ 



exp 



i { 2m + T 



xir 
M 



(53) 



where 



r 



n 



n 



0,1, 



'N + l 



N -1 



.N 
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and 



N 



(2y+l)^. 



n = 0, 



n = l,2,---, 



N -1 



N 



N 

n = — , for even only, 



(2y + l)^ 



n 



N 
~2 



Substituting these expressions in (fTT|) . separating out the summation for 
n = 0, and making use of the identity 



sm 



{2y+l)(- + n) 



N 



-lycofi {2y+l) 



nir-i 



we obtain after some reduction 

M-l N-1 

-R{MxAr}(ri, = X! X! 

m=0 n=0(m,n)7^(0,0) 



I.Klein „/.Klein 

{m,n);{xi,yi) ^ {m,n)\(x2,y2) 



r 



M-l Af-1 



M 



+ E E 



'o n=l ^(m.,n)(Tn) 



Klein 



Klein _ 

(m,n);{xi,yi) 1^ (m,n);{x2,y2) 



r 
N 



Xi - X2 



{Xi - X2f 



M 



A 



N 



riV-li 

1 M-l i^-l 



+ ^1 +^2 -2(-l)( ^^-^^)"CiC2 cos [ 2(a;i - X2)e^(r) ] 



MN 



T=0 m=0 n=l 



^{m,n) 



(54) 



where 
An 



T\ TT 

m H — — 
2/M 



M-l 

E 



1 _ (_l)?/i-?/2 cos[2(xi-X2)e„(l)] 



0, AT = odd, 



A(m,Ar/2)(l) 



iV = even 



(55) 
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and Ci = cos[{yi + l/2)mT/N],i = 1,2, as defined in (pUj) . 

Example 10: The resistance between nodes (0, 0) and (3, 3) on a 5 x 4 
(A^ = even) Klein bottle with r = s is computed from (j5l|l as 

figs, ({0.0}, {3, 3}) . (iL^I^H). 

= (0.654149...) r, (56) 

where the three terms in the first line are from the evaluation of corresponding 
terms in (j54j) . The result is to be compared to the corresponding value for 
the same 5x4 network under the Mobius boundary condition considered in 
Example 9, which is the Klein bottle without periodic boundary connections. 



9 Higher-dimensional lattices 

The two-point resistance can be computed using (fTT|) for lattices in any 
spatial dimensionality under various boundary conditions. To illustrate, we 
give the result for an M x x L cubic lattice with free boundary conditions. 

Number the nodes by {m, n.i}, 0<m<M-l,0<n<N-l,0<i< 
L — 1, and let the resistances along the principal axes be, respectively, r, s, 
and t. The Laplacian then assumes the form 

where Tf^'^'^ is given by (jl9|) . The Laplacian can be diagonalized in the three 
subspaces separately, yielding eigenvalues 

\m,n/) = 2r"^(l - cos6'^) + 2s"^(l - cos0„) + 2t"^(l - cosa^), (57) 

and eigenvectors 

where il^^J is given by ipUl) and = iir/L. It then follows from (fTT|) that 
the resistance -Rfree between two nodes ri = {xi,yi,Zi) and r2 = {x2, 1/2,^2) 
is 

A/-1 N-lL-l 

-R{MxArxL}(ri,r2) = X! X! 51 ^{m,n/) 

m=0 n=0£=0(m,n/)^(0,0,0) 

If f ^ 

' (58) 
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The summation can be broken down as 

r,free / \_ ST^ ST^ \^{m,n,e);{xi,yi,zi) r(m,n,£); (0:2, 1/2,^2) 

^{MxNxL}[''^li^2) — 2^ 2^ 2^ 



m=l n=l £=1 ^{m,n/) 
1 



1 1 

+YR{MxN}{{xuyi}Ax2,y2}) + jjR{NxL}{{yu zi}Ay2, ^2}) 



-^R{LxiM,x2) - ^Rl\7.,}iyi,y2) - -^R^nxiM^z^). (59) 

All terms in (j^^ have previously been computed except the summation in 
the first line. 

Example 11: The resistance between the nodes (0,0,0) and (3,3,3) in 
a 5 X 5 X 4 lattice with free boundaries and r = s = t is computed from 
as 

R'- no 01-13 3 31) - r 327687658482872 \ 
%x5x4}(|0,0,0|, |3, 3, 3|) - 352468567489225 

= (0.929693...) r. (60) 



Example 12: The resistance between two interior nodes ri and T2 can 
be worked out as in Example 4. The result is 

2 /'27r /'27r /'27r 

-Roo(ri,r2) = ,^ ,„ / d(j) I d9 da 



X 



(27r)3., 

1 — cos (^Xi — X2)0 cos(?/i — 2/2)0 cos(zi — 2:2)0; 



r ^(1 — cos6') + s ^(1 — cos0)+t ^(1 — cosa)/' 
which is the known result jl]. 



10 Summation and product identities 

The reduction of the two-point resistances for one-dimensional lattices to 
the simple and familiar expressions of and (jH^ is facilitated by the 
use of the summation identities ()27|) and ljH^ . In this section we extend the 
consideration and generalize these identities which can be used to reduce the 
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computational labor for lattice sums as well as analyze large-size expansions 
in two-and-higher dimensions. 

We state two new lattice sum identities as a Proposition. 

Proposition: Define 

'V„t'o<:oslii-cos(a!f)' ' 

Then the following identities hold for I > 0, N = 1,2, ■ ■ ■ , 



^ cosHN-i)l 1 
^ (sinh/)sinh(M) 



1 , l-(-l)^ 



+ 



sinh^/ 4cosh2(//2) 



< £ < 2iV, 

(61) 

- (smht£M/2) - «^^<^- 

Remarks: 

1. It is clear that without the loss of generality we can restrict £ to the 
ranges indicated. 

2. For £ = and / 0, /i(0) leads to ^ and hiO) leads to (j^. 

3. In the N oo limit both and (jU^ become the integral 

cos(W) ^^^^ j^3, 



TT JO cosh Z — COS 9 sinh / 

4. Set £ = in (jF)T|l . multiply by sinhZ and integrate over I, we obtain 
the product identity 



Yl I cosh / - cos — J = (sinh Nl) tanh(Z /2) . (64) 

n=0 ^ ^ ^ 



5. Set £ = in ()62|1 . multiply by sinhZ and integrate over /. We obtain 
the product identity 

^-1 / 2n7r\ 

[] (^cosh/ -cos^j = sinh^(M/2). (65) 
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Proof of the proposition: 

It is convenient to introduce the notation 

^ iV l + a2-2acos^„' 



a < 1, a = 1, 2 



n=0 



so that 



I^{£) =2aSaii), a = e-K 
It is readily seen that we have the identity 

1 r 



SJl] 



2a 



(l + a")5„(0)-l 



(66) 



(67) 



(68) 



1. Proof of dnH): 

First we evaluate 5'i(0) by carrying out the following summation, where 
TZe denotes the real part, in two different ways. First we have 

Af-l 



1 1 



1 



a e' 



AT ^-^ 



N n=o |1 - ae'^^ 

1 1 — a cos 9n 

N ^0 l + a2 -2acos^„ 
5i(0)-a5i(l) 

l\l + {l-a^)Si{0) . 



(69) 



Secondly by expanding the summand we have 



1 1 

n=0 



N-l oo 



a e 



a e 



n=0 ^=0 



and carry out the summation over n for fixed £. It is clear that all C = 

even terms vanish except those with i = 2mN,m = 0,1,2,... which yield 

Eoo 2mN 
m=0 " 



1/(1 - a^^). For ^ = odd = 2m + 1, m = 0,1,2,... we have 



Af-l 



n=0 ^ 



l-(-l) 



2m+l 



oi{2m+l)iT/N 
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after making use of (j25|) . So the summation over I = odd terms yields 
N-^ E^=o a^™'*'^ = - «^), and we have 



Af-l 



7^e V = = - 

^nl-a e'^- 1 - N{1 - a^) 



+ 



(70) 



n=Q 



Equating (j^Hl) with (fTUI) we obtain 



^i(O) 



l-a2 



1 + a 



+ 



2a 



(71) 



To evaluate 5*1 (£) for general £, we consider the summation 



1 ^'l-(ae^'"Y 



TZe-y — 
iVfr^n 1 



n=0 



a e 



AT ^ 



n=0 



l_a^e^^^")(l-ae 
1 1 — a e*^" I ^ 



5i(0) - a5i(l) - a^Si{£) + a^~^^Si{£ - 1)(72) 



where the second line is obtained by writing out the real part of the summand 
as in ()69|). On the other hand, by expanding the summand we have 



1 ^1 1 - (ae'^- 



N-l £-1 



TZe-Y — 
N f4 1 



n=0 



a e 



TZei- V V 



m inmn/N 



a e 



n=0 m=0 

1 e-i 



m=l 



1 



1 + 
1 + 



a(l — a ) 

ail - a^-n 



i = even < 2N 



odd < 2N, (73) 



Ar(l -a2) ' 

where again we have used (j2SI)- 

Equating (fTSj) with (f?^ and using (jUHj) and (f7T|). we obtain the recursion 



relation 



where 



5^(£) - a 5jv(^ - 1) = A + 5^ 



,2 AT 



1 — a 



2N ' 



5, 



N{l-a^) ■ 



(74) 



(75) 
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The recursion relation (j74j) can be solved by standard means. Define the 
generating function 



a = 1,2. 



(76) 



Multiply (j74j) by and sum over i. We obtain 



\-a-H iV(l -a2)(l -t2)- 



(77) 



This leads to 
Gi(t) 



1 - at 



SM + - 77 + 



t + at^ 



l-a-H iV(l-a2)(l-t2) 



1 



(1 -a2)(l-a2^) 
1 



1 a 

+ 



2N 



1 — at 1 — a ^t 
1 



2A^(1 -a)2(l -t) 2A^(l + a)2(l + t)^ 



from which one obtains 



^i(^) 



a + a 



2N-£ 



+ 



(1 -a2)(l-a2^) 2A^(l-a)2 2N{l + ay 



a + a 



2N-i 



+ 



4a 



1 



(l-a2)2 (l + a2) 



(78) 



(1 -a2)(l -a2^) 2iV 

It follows that using Ii{£) = 2aSi{£) we obtain (jHTjl after setting a = e~K 
QED • 

2. Proof of dni: 

Again, we first evaluate 5*2(0) by carrying out the summation 

N-l -| 

(79) 



1 1 



n.=0 



a e 



a < 1 



in two different ways. First as in (j69|) we have 

N-l 



1 '"-^ 1 



a e 



2 
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l + (l-a')52(0) 



(80) 



where 5*2 (i) is defined in . Secondly by expanding the summand we have 

1 N-l 1 1 N-1 oo -1 

^ E rrips: = ^ E E = (81) 

where by carrying out the summation over n for fixed i all terms in ()7Up 
vanish except those with i = mN, m = 0, 1, 2, ... Equating (jHT|) with ()80|) we 
obtain 



and from (|UHj) 



We consider next the summation 

[ae ) 

iV ^0 1 - a e*2e„ ^ ^ 

Evaluating the real part of the summand directly as in (ff^ . we obtain 

ne-Y: , 2 .2 J = ^2(0) - aS,{l) - a'S,{i) + a'+'S,{i - 1). (84) 
Secondly, expanding the summand in (jH!?jl we obtain 



-1 gi2m7r 



m i2TTmn/N 



1 

iV 



m=l 



^i2rmT/N 

1 m < £ < A^. (85) 



Equating and (|S^ and making use of (jH2I) for 5*2(0), we obtain 

1 — a 



S2{i) - aS2{i - I) = (86) 
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The recursion relation (j86|) can be solved as in the above. Define the 
generating function G2(t) by (f7B|) . We find 



G2it) 



1 



^2(0) + 



1 — at 



(l-a^)(l-a-it) 



1 



(87) 



{1-a 



2^) [l - at 1-a-H 



from which one reads off 



f I N—f 

a + a 



(88) 



(l-a2)(l-a- 



2N 



Using the relation /2(^) = 2aS2{i) with a = e we obtain (jU^ . QED • 
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Figure Captions 

Fig. 1. A network of 4 nodes. 

Fig. 2. A one-dimensional network of M nodes with free ends. 

Fig. 3. A one- dimensional network of M nodes with periodic boundary 
conditions 

Fig. 4. A 5 X 4 rectangular network. 

Fig. 5 (a) The schematic plot of an Mobius strip, (b) The schematic plot 
of a Klein bottle. 
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Fig.l 
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